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^^ , We introduce and analyze a general one-dimensional model for the description of 

^^ \ transient patterns which occur in the evolution between two spatially homogeneous 

^^ \ states. This phenomenon occurs, for example, during the Freedericksz transition in 

nematic liquid crystals. The dynamics leads to the emergence of finite domains which 

are locally periodic and independent of each other. This picture is substantiated by 

J^ , a finite-size scaling law for the structure factor. The mechanism of evolution towards 

the final homogeneous state is by local roll destruction and associated reduction of 

local wavenumber. The scaling law breaks down for systems of size comparable to 

the size of the locally periodic domains. For systems of this size or smaller, an 

apparent nonlinear selection of a global wavelength holds, giving rise to long lived 

periodic configurations which do not occur for large systems. We also make explicit 

J2 \ the unsuitability of a description of transient pattern dynamics in terms of a few 

Fourier mode amplitudes, even for small systems with a few linearly unstable modes. 
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I. INTRODUCTION 

The most simple situation considered in the context of pattern formation studies [jl| is 
the one in which an homogeneous stable steady state of a system becomes unstable at a 
threshold value of a control parameter, so that beyond threshold the new stable state is 
time independent and with a well defined spatial periodicity. In such situation it is first 
generally aimed to describe some static properties, such as the threshold value of the control 
parameter, possible wavelengths of the pattern and possible higher-order bifurcations. Other 
set of interesting questions is associated with the transient dynamics of the pattern formation 
process given an initial unstable homogeneous state 0. A different physical situation is the 
one of transient pattern formation during the temporal evolution between two homogeneous 
steady states: In this case an homogeneous stable steady state becomes unstable beyond 
a given threshold and the transient evolution to the final homogeneous stable steady state 
involves a process of pattern growth and decay. Most of the usual mathematical techniques 
used to describe the first situation seem to fail for problems involving transient pattern 
dynamics. 

Transient pattern formation is well documented from the experimental point of view 
for different instabilities in nematic liquid crystals P,^. A typical situation is the mag- 
netic Freedericksz transition in which for a large enough applied magnetic field the nematic 
director does not reorientate homogeneously but a striped pattern with a characteristic 
wavelength emerges. This pattern can last from minutes to hours depending on the specific 
geometry and the material, but it finally disappears leading to the homogeneously-reoriented 
final equilibrium state. For most geometries of the system, the pattern admits a good de- 
scription by one-dimensional models. A good summary of different situations considered in 
this context is given in Ref . . Generally speaking the transient pattern is associated with 
the coupling of the director field with hydrodynamic variables so that the fastest response 
to the applied field is not homogeneous in space . A mechanism of wavelength selection 
based on this idea of fastest response has been proposed ^-0]: the well-defined observed 
periodicity has been associated with the mode of fastest growth and its dependence with 
the applied magnetic field has been considered theoretically and experimentally. However, 
a nonlinear mechanism of wavelength selection has also been invoked and substantiated by 
some experimental and restricted numerical studies |§| . The process of pattern formation 
in nematics can be described in detail through the full set of nematodynamic equations 
which have been consistently formulated also in the presence of thermal fluctuations 
0. Such equations have been discussed in a variety of situations 0,^J^,|]. These detailed 
discussions of a rather complicated set of equations might hide some general features of 
the problem of transient pattern formation and development. Our aim in this paper is to 
propose and analyze a generic model for one- dimensional transient pattern formation which 
describes general features of these problems. It is obvious that a precise comparison with 
experiment would require the consideration of some specific details. Nevertheless, we hope 
that general aspects as wavelength selection (if such selection does occur), mechanisms of 
pattern evolution and the important issue of finite size effects can be understood from this 
general model. 

The model we analyze is defined by the following equation for a scalar variable 6{x,t): 

e{x,t) = {a-dj){dje + c0-be^) . (i) 



The dot denotes temporal derivative. The hnearized version of this equation in Fourier space 

e, = uj{q)e, (2) 

involves an amplifying factor 

uj{q) = {a + q^){c-q^) . (3) 

As seen in Fig. 1, this factor is such that for c > there is a range of unstable modes 
for which uj{q) > 0. If a < the range of unstable modes does not include the mode 
g = and the linear regime is qualitatively the same as for the well known Swift-Hohenberg 



equation |10| used to describe the formation of stationary patterns. We are here interested 
in the situation in which a > for which the range of unstable modes is {—qcQc), with 
qc = c^/^. If in addition, a < c, the mode of fastest growth becomes different from zero: 
Qm = ((c — a)/2)^/^. This implies the existence of an instability that for < a < c involves 
the linear growth of a pattern with a characteristic periodicity given by g^- Throughout the 
paper we fix c = 1 and study (|ID with periodic boundary conditions. 

Equation (|I|) can be motivated as an approximation to the complete nematodynamic 
equations describing the magnetic Freedericksz transition in a nematic. Under ordinary 
assumptions for the twist geometry one obtains coupled equations for the angle orientation 
6{x) of the director and a component of the velocity field, both in a plane in the middle of 
the sample. The approximation of negligible inertia permits the elimination of the velocity 
variable. This gives rise to an effective wavenumber-dependent viscosity which appears as a 
g-dependent kinetic coefficient Q. For small deformations 6 and in the small wavenumber 
limit [|n|,|T^ one recovers equation (0), where the factor (a — dj^) is the remaining part of 



the effective viscosity [[T3 |. 



Within the context of model (Q), we address in this paper three general questions associ- 
ated with transient patterns dynamics in the intermediate nonlinear regime after the initial 
pattern emergence and before the late stages in which it disappears. The first question is the 
domain of validity of the linear theory and the possible existence of a wavelength-selection 
principle in the nonlinear regime. The second question is the characterization of the mech- 
anism governing pattern evolution. Thirdly we examine finite-size effects which could be 
preponderant in the question of a nonlinear selection principle. Our results indicate that 
for large systems a wavelength is initially selected in the linear regime but, on the average, 
it then changes monotonically in time approaching a final homogeneous state. During this 
evolution the system can be described as composed of several regions which evolve indepen- 
dently of each other. A local mechanism of roll destruction operates in such regime. This 
description is substantiated by a scaling law for the structure factor of systems of differ- 
ent size which reveals the presence of uncorrelated regions of a characteristic length. This 
length is essentially time-independent in the regime of dynamical evolution examined. For 
systems of size comparable or smaller than this length an apparent nonlinear selection of a 
wavelength occurs: configurations of well defined periodicity, which is not the one linearly 
selected, last for very long times, before evolving to the final homogeneous state. The reason 
why such periodic configurations are long lived is that they correspond to unstable station- 
ary solutions which, for small system sizes, are approached in the initial transient regime. 
For large system sizes such configurations are not approached in the dynamical evolution 
from a typical initial configuration. A characteristic of small systems is that a very small 



number of modes are linearly unstable at t = so that a description in terms of coupled 
ordinary differential equations for the amplitude of a few modes seems natural. We show 
that such description can give qualitatively wrong results. The basic simple reason for this 
fact is that such modes are not stabilized by nonlinear terms in the full equation, because 
they are not associated with stable solutions. 

We finally wish to mention the relation of our analysis with another well-known physical 
situation in which transient patterns occur, namely spinodal decomposition WM. The pro- 



cess of phase separation of a binary mixture quenched to a temperature below the critical 
temperature also displays a transient pattern with a time dependent characteristic length. 
This analogy was pointed out many times 0. However, in addition of differences in time 
scales P] there is another important difference between this problem and the one of transient 
patterns in nematics and it is that the dynamics of spinodal decomposition is constrained 
by a conservation law. The fact that a 7^ in (|ID indicates that there is no conservation 
law for the spatial integral of 6{x,t). One-dimensional {d = 1) spinodal decomposition is 
physically artificial and consequently has not been studied in great detail from the physical 
point of view. The initial linear regime (Cahn-Hilliard-Cook theory WM) and the very late 



stages of phase separation ||T^ are well understood, but no detailed study seems available for 



the d = 1 intermediate nonlinear regime in which we are here interested. The initial PJTB[ 
and very late stages [0] of pattern dynamics in the magnetic Freedericksz transition have 
been already studied in some detail by analogy with the problem of spinodal decomposition. 
Standard nonlinear theories of spinodal decomposition such as that of Langer, Bar-on, and 
Miller |l^ do not seem applicable in d = 1: since they are based on the competition be- 



tween locally ordered equivalent stable states they do not include the competition between 
different wavenumbers |1^ which seems essential in d = 1. The other classical problem of 



the dynamics of phase transitions, namely the dynamics of an order- disorder transition |T^ 
does not involve a conservation law but it neither has a periodicity selected by the dynamics 
of the process. 

The outline of this paper is as follows. In section 2, we discuss the main characteristics 
of the model given by (0) and its stationary configurations, then we compare it with other 
related models. Section 3 describes our numerical results for large systems. Finally the 
description of a small system in terms of the amplitudes of a few modes is discussed in 
section 4. Throughout the paper we restrict ourselves to d = 1 and we neglect thermal 
fluctuations. The role of fluctuations and two dimensional effects will be analyzed elsewhere. 

II. A MODEL FOR GENERIC ASPECTS OF TRANSIENT PATTERN 

FORMATION 

A first important property of the model given by equation (|l|) is that it can be written 
in a potential form 

where the kinetic coefficient F is the operator 

T = a-d^\ (5) 



and F[e] is 



m 



dx 



\e{xf + ^-e{xf + \\d.e{x)\' 



(6) 



From ([l|) we can show that F is a good Lyapunov functional, in the sense that dF/dt < 

dF[e] 



dt 
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The last inequahty holds because F is a positive definite operator (if a > 0), as can be seen 
from its expression in Fourier space. 

The picture of the evolution is then that the system evolves, continuously decreasing 
F[6], until a minimum of F[6] is found, thus stopping the evolution. Note that such minima, 
stationary solutions of {^, are independent of F and of the parameter a and they are solutions 
of the simpler equation 



dje{x) + 9{x) - he{xf = 



(8) 



An independent demonstration of the fact that with periodic boundary conditions the only 
stationary solutions of (|l|) are those of (H) can be set up by writing (|I]) with 9{x) = as the 
set of two equations: 



d^'^9{x) + 9{x) - h9{xf = y{x) 



(9) 



Ty{x) = 



(10) 



Eq. (10) is a linear second order ordinary differential equation whose general solution is a 
linear combination of two exponentials. Only the combination leading to y{x) = satisfies 
periodic boundary conditions, so that (^ reduces to (|). The qualitative features of all the 
stationary solutions of (0) can be discussed by writing (||) as [119 



d^9{x) 



dV{9) 



(11) 



dx"^ d9 

which resembles a Newton equation for the motion of a particle of unit mass in a potential 



V{9) 



1^2 _ ^^4 

2 4 



(12) 



the role of 'time' being played by the coordinate x. From this analogy, the bounded solutions 
of (^ can be classified in three types: 

a) The uniform solutions 9{x) = l/\/b and 9{x) = —l/\/b. 

b) The uniform solution 9{x) = 0. 



c) A family of solutions represented by nonlinear oscillations in the potential V{6). The 
maximum possible frequency corresponds to oscillations of small amplitude around 
^ = 0, being its period 27r. The minimum frequency corresponds to trajectories in 
which 9{x) remains mostly near ±l/v6, with short excursions (domain walls) linking 
both states. In summary: there are periodic solutions to (^, which we denote by 
ipq{x), each one containing a different fundamental wavenumber q and its harmonics, 
with < g < 1. 

An important question is the linear stability of such stationary solutions. In order to 
consider this question, Eq.(^ does not contain enough information and the full dynamical 
equation (|l|), linearized around the stationary solution being checked, is needed. Lineariza- 
tion around the uniform solutions is immediate and it is found that the solution 6' = is 
linearly unstable, and both 9 = 1/Vb and 9 = —1/Vb are linearly stable. These are also 
the absolute minima of the functional F[9], so that they represent the stable equilibrium 
phases. The analysis of the stationary periodic solution ipq{x)^ of fundamental wavenumber 
q, is performed with the introduction of 9{x,t) = ipq{x) + A(x, t) and linearization in A. 
The resulting equation is 

A(x, t) = r [l - ?,hil)q{xf + d/\ A(x, t) . (13) 

The general analysis of this linear equation with periodic coefficients requires of Bloch or 
Floquet theory |^. A simplified situation was considered in |]1T| for the case in which q 



is small, so that the solution consisted basically of domains of the stable phases separated 
by thin domain walls. In that case it was found that the periodic solutions were linearly 
unstable. It can be generally shown that all the periodic solutions are unstable by studying 
its stability with respect to a uniform perturbation. The argument is as follows: Let us 
introduce the initial perturbation A(x,t = 0) = Aq, Vx, and consider the initial time 
t = 0"*", when A{x,t) ^ Aq. Let be a; = one of the places in which ipg = 0. Near x = 0, 
ipq{xY will have a positive parabolic shape, so that d'^[il'q{x ~ 0)^] > 0. Introducing this in 
Eq.(p!3l) we find that sign[A(x ~ 0, t = 0"*")] = sign[Ao], showing the instability of ipq because 
a uniform initial perturbation grows. The consequence of the instability of all the periodic 
solutions is that none of them can represent the final state of the evolution, as long as a 
non-zero amplitude is given in the initial condition to the mode with wavenumber g = 0, or 
if noise is present in the system. 

After this summary of the general properties of Eq. (^ it is interesting to compare them 
with the properties of other related models studied in the literature. To this end we write 
(0) (with c = 1) as 

9 = -dj9 + (a - l)dj9 + a9- ab9^ + bdj9^ . (14) 

For < a < 1, the uniform solution 6* = is unstable and the linear analysis predicts 
the growth of modes with wavenumber q y^ 0. It is then natural to relate this equation with 
others for which a periodic pattern grows from the unstable uniform solution. An archetypal 



example of such equations is the Swift-Hohenberg equation |]T0 

e=\j^-{l + d,^f] 9 - b9'' . (15) 



The linear stability analysis of this equation leads in fact to the same linear growth spectrum 
than Eq.(|l|) except for an important difference in sign in the regions around g ~ 0: the 
modes in this region slowly grow in our model (a > 0) whereas they are damped in the 
Swift-Hohenberg case (due to the fact that 7^ is positive). Other aspects of the initial stages 
in pattern development are qualitatively similar for both models. Another more fundamental 
difference is that the Swift-Hohenberg equation admits a family of stable periodic solutions, 
one of which will give the final state of the evolution, whereas (|ip admits no other stable 
solutions than the uniform 9 = ±l/^/b. This difference comes from a combination of the 
different sign of 7^ versus —a and of the additional term dx'^9{x)^ in (|T1|). Thus, the evolution 
at late times will be completely different in both cases. 

These differences have important methodological consequences: in order to study pattern 
formation in cases exemplified by the Swift-Hohenberg equation, a common strategy is to 
take as a small parameter the range of unstable modes around the most unstable one, which 
is small near a bifurcation point, and then obtain a nonlinear equation for the amplitude of 
the most unstable mode. The form of this equation is greatly determined by the symmetries 
of the problem, and by the assumption of being the first step in a uniform expansion. 
This strategy can not be applied to our problem, because the characteristic shape of the 
linear instability spectrum in (^ with a fastest growing mode g^ 7^ is only obtained for 
c > a which is far enough from the bifurcation point (c = 0), and the band of unstable 
wavenumbers is as large as the wavenumber of the most unstable mode (because the mode 
with g = has to be included in the description) . Then bifurcation theory and normal forms 
are of no much help in our problem. In addition, the fact that the mode of fastest growth 
is not associated with a stable solution precludes the use of approximations based in the 



saturation of the linearly fastest growing mode, such as those in |^ . 

Another class of models with which it is natural to compare our model is the one repre- 
sented by the Fisher-Kolmogorov equation, also known as Ginzburg-Landau equation for a 



real variable, or, with an added noise term, model A of critical dynamics |T^,|22|: 

e{x,t) = dje{x) + e{x)-be{xf . (i6) 

The stationary solutions of this equation are exactly the same as in our model, and the 
only stable solutions are, as in our case, the uniform 9 = ±l/y/b. In fact, the analysis in 
|TT| shows that at very long times, the dynamics of the domain walls in (|l|) is the same 



as in (|TC|). The main difference is, however in the conditions created by the initial linear 
instability if a < 1. The fastest growing mode in (0) is the one with wavenumber q = 0, 
which corresponds also to the final state. Then theories such as those in PT| are good 
descriptions of the evolution for all times. 

In some sense, the time evolution of our model is a crossover between a linear behavior 
close to that in the Swift-Hohenberg model, and final stages similar to those in (|16|). Perhaps, 
this is why the closestly related model is the one represented by the Cahn-Hilliard equation, 
describing spinodal decomposition in binary mixtures and alloys. It is also known, when 
a noise term is added, as model B of critical dynamics 0,^. Formally, this model is 
obtained by putting a = in Eq.(|I|). In this case we have also an initially periodic structure 
which coarsens in time to approach q = 0. The main difference with our model is that the 
spatial integral of 6{x) is conserved by the Cahn-Hilliard equation, so that a completely 
uniform solution can not be approached unless J dx 6{x,t = 0) = 0. In the generic case. 



the final state is the coexistence of two domains of the stable phases and not only one as 
in our case. Thus, the final stages of evolution should be very different in both models 
pn| . The fact that, in addition to the fastest growing mode, the mode with q = is also 



linearly unstable in @ implies a wide range of unstable wavenumbers in the initial regime, 
leading to a wide spectrum during the nonlinear stages. Time- dependent configurations 
do not approach closely to any of the unstable periodic solutions. The consequence is that 
theories such as that of Langer [^ assuming that the system is close to one of the stationary 
unstable solutions at each time, will be only of certain usefulness at extremely long times 
n|, where the mode with q = will be the dominant one. 



III. NUMERICAL STUDY OF TRANSIENT DYNAMICS FOR LARGE SYSTEMS 

The time evolution of 6{x, t) from an initial condition close to the unstable steady state 
6{x) = has been calculated by solving numerically equation (|l]) on a grid of N points. In 
the remaining part of the paper, we fix the value of the parameters in ([l|) to a = 0.002, 
6 = 3 and c = 1 as appropriate for typical values of the parameters in the nematodynamic 
equations ||,|l^. In this case g^ = 0.7. We have used a centered finite-difference scheme up 
to order (rfx)^ to approximate the spatial derivatives. A predictor-corrector method with 
one step has been used to determine 9a.tt + dt. A suitable value for dx has been determined 
integrating the equation in the linear regime and comparing the growth rate of the unstable 
modes obtained numerically with the one calculated analytically. A value of dx = 0.25 has 
been thereby chosen. The most unstable mode has a wavelength of Xm = 27r/gm ~ 8.98 
which corresponds to approximately 36 grid points. The length of the system is L = Ndx. 
We have considered a range of system sizes from L = 64 to L = 256 and we have always 
taken periodic boundary conditions. The time step used is 10~^. For dt larger than 0.004 
numerical instabilities were observed. For values of dt ranging from 2 ■ 10~^ to 5 ■ 10"^ the 
discrepancies in 6{x) after 5 units of time of integration were smaller than 10~^. 

The initial condition is written as 

^{x,0) = - + e^sin{qnX + (j)qJ (17) 

where the sum is over modes g„ = 2iTn/L and it has been usually taken to run only over the 
unstable modes g„ < qc] 0g„ is a random phase and e is a small amplitude assumed equal 
for all of the modes and whose value was arbitrarily set equal to 2 ■ 10~^. To obtain different 
initial configurations of 6{x) the set of random phase shifts (pq^ was changed but not the 
amplitude e. Hence, when talking about a different initial condition we mean a different set 
of random phase shifts. 

In addition to our transient dynamics study, we have also examined the existence of 
periodic stationary solutions (type 'c' in the previous section). Starting with a configuration 
of the form 9{x) = e sin{qnx) involving a single mode g„ < qc with an amplitude e = 
2 ■ 10~^, the pattern develops with the growth of its harmonics until a stationary pattern 
is obtained. We know that this pattern is unstable and the velocity of its decay has been 
tested numerically by adding a small amplitude to all the modes with q smaller than qc- 
The decay was found to be always extraordinarily slow. This means that these unstable 
stationary solutions can be long lived. When the initial condition (|17|) is used, the pattern 
develops in a way that none of these periodic stationary patterns is approached during the 
transient dynamics provided the system size is large enough (see next section). 

Our results for the transient dynamics study are summarized in figs. 2, 3 and 4 for the 
evolution of the configuration 6{x, t), the associated structure factor and the number of rolls 
of the pattern respectively. General features of the time evolution which manifest themselves 
in specific ways in these figures are the following: A linear and a nonlinear regime of evolution 
can be clearly identified. In the linear regime the pattern is formed. In the nonlinear regime 
no mechanism of wavelength selection occurs, but domains of a characteristic size exist. 
These domains include several rolls and evolve in time in a way essentially independent of 
each other. 



In Fig. 2a, the time evolution of the pattern 6{x) for a particular initial condition is 
shown. The periodicity related to the linearly most unstable mode becomes apparent on 
6{x) in the initial stages. Afterwards, we observe that rolls disappear continuously. The 
presence of regions of different periodicities can be observed. For example, at t = 40 around 
7 of these regions are distinguished. A systematic method of identifying such regions is to 
find the power spectra of small regions in the pattern, and identify the maxima in these 
local spectra with a dominant local wavenumber 0. Explicitly, we have multiplied the 
configurations in Fig. 2a times a Gaussian of width a = 10, unit height, and centered at 
X. Then we have calculated the maxima in the power spectra of such localized patches 
as a function of x. Figure 2b shows the result of one of such analysis, at time t = 40, 
identifying an x-dependent dominant wavenumber. At this time the average wavenumber 
has already deviated from the linearly most unstable one {qm = 0.70) so that one has entered 
the nonlinear regime. In Fig. 2a we can also see that, at each time, a roll is disappearing 
in that region whose local dominant wavenumber is largest. Note that the disappearance 
of a roll occurs locally since this does not affect other regions (compare for example 6{x) 
at t = 100, 200 and 300). At the longest times arrived, there are regions were 6 is already 
close to the value of the uniform solution 9 = ±6^2 ^ 0.58. After this time, it is expected 
that the pattern evolves towards configurations where regions of ^ = 6~^'^ and 9 = —h~^''^ 
coexist separated by walls. The separation between these walls will be rather large and the 



pattern evolution should be described by the theory in |TlI . 

The structure factor S{q,t) associated with 9{x,t) is defined from the discretized config- 
uration {9{xn = ndx,t), n = 1, ..., A^} as 



S{l,t) = ^ 



Y,e"--e(x^,t) 



n 



(18) 



The vertical bars denote the modulus of a complex number. Note that a normalization 
factor 1/A^ has been included. With this choice of normalization S{q) is independent of 
system size for a uniform configuration ^(x„, t). It is also independent of N during the linear 
regime. The allowed values of q are of the form n dq with n an integer between —N/2 
and A^/2 and dq = 2tt/ L. Figure 3a shows the time evolution of S{q) averaged over 50 
independent initial conditions of the form (plT]) . Figure 3b shows the time evolution of the 
area A[t) = J dqS{q,t) for 4 of these independent initial conditions. Since linear theory 
predicts that the structure factor is independent of the initial phases (pg^ , the time at which 
the different curves in Fig. 3b begin to separate signals the end of the linear regime. This 
happens around t ~ 20. During the linear regime the structure factor in Fig. 3a is shown 
to grow with a maximum around the linearly most unstable mode qm = 0.7. 

After the linear regime, and reflecting a continuous elimination of rolls, the maximum 
of S{q,t) shifts towards small g's as times goes on (see Fig. 3a). This continuous drift 
characterizes the elimination of rolls and it eliminates the idea of a nonlinearly selected 
wavelength. During the nonlinear regime the structure factor develops important contri- 
butions for short and long wavelengths which indicate the existence of strong competition 
between many modes. The existence of different domains in the configuration 9{x, t) can be 
characterized from the structure factor by defining a correlation length l^ = 2it/w, where 
w is one half of the width of the peak of S{q,t) at its mean height. This length has been 
plotted in Fig. 3c for a system of size L = 128. It is around 1/3 of the system length for this 
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system size. Hence, we talk about approximately 3 uncorrelated zones in the system which 
evolve independently. The validity of this view is enhanced by the evolution observed in 
systems with L = 256, as discussed above, and also with L = 64: w turns out to be roughly 
independent of system size and time (during the nonlinear regime we consider here). The 
domain size is probably determined by the interplay between the intensity e of the initial 
condition and the shape of uj{q). 

Figure 4 shows the average number of zeros per unit length of 6{x) as a function of time. 
This quantity identifies the number of rolls per unit length of the pattern and it measures 
the 'average' wavenumber in the system. The inset in this figure represents the average over 
initial conditions of the number of zeros per unit length vs a mean wavenumber defined as 
(g) = J qS{q,t)dq/A(t). The upper part of the curve in the inset corresponds to the very 
early initial regime. Beyond this regime the number of zeros is linearly related to (q) so 
that a description in terms of any of these quantities is equivalent. The number of zeros 
per unit length is seen to grow during the linear regime of pattern emergence reaching a 
value around 0.22 which corresponds to the global wavenumber selected by fastest linear 
growth, qm = 0.7. This number of zeros remains constant for a little while beyond the end 
of the linear regime at t ~ 20. Later the average number of zeros decreases monotonically 
in time making again clear that we can not identify any nonlinearly selected wavenumber. 
This result is different to the one found in Ref. 8 for a similar system. We will show in 
the next section that an apparent wavelength selection might occur due to finite size effects. 
We have attempted to obtain a growth law for the decay process in Fig. 4. The decay is 
certainly slower than the one given by a power law. A relation (q) = 21.3 — 1.27 In t seems 
to fit the curve for t > 90. Nevertheless, a decay (Int)^ could also be fitted. Determining 
the exact law and an accurate exponent will need of a considerable increase in the statistics 
and it is not our goal here. We note that a logarithmic decay is the one expected from 
one-dimensional domain- wall dynamics ||ll],|l5l, but the largest times in our simulations are 



still far from the regime in which the pattern is composed of domains of the equilibrium 
phases separated by thin walls. 

An important question in relation with the existence of domains of a characteristic size 
is the dependence of the transient dynamics on the system size L. It is already seen in Fig. 
4 that the evolution of the number of zeros per unit length is essentially the same for two 
system sizes. We have checked that the time scales of evolution are independent of system 
size for L > 50. Taking advantage of this fact, a systematic study of the dependence on L 
of the temporal evolution of the structure factor can be carried out by the analysis of the 
evolution of the area A of S{q,t). This area is shown in Fig. 5a for systems of different 
size. From t = 75 to the longest time considered {t = 200), a relationship of the form 
A{L,t) = L°'f{t) is satisfied, with a = —0.49 ± 0.01 ~ 1/2, for system sizes large enough. 
Figure 5b shows the validity of this scaling behavior. The combination of this result with 
the existence of size-independent time-scales of evolution implies that the dependence of the 
structure factor on system size factorizes out: 

S{q,t,L)=L-"^F{q,t) (19) 

It should be stressed that this factorization is not a trivial consequence of the normalization 
factor in (p!8|). This finite-size scaling form contrasts with the one found in the dynamics 
of order-disorder transitions and of spinodal decomposition p^ in d > 2, where the scal- 
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ing function F depends on system size as F{qL,t/L^), and the exponent a is instead of 



— 1/2 1 24]. The reason for these differences can be elucidated by noting that the exponent 
a = —1/2 is a manifestation of the fact that the system is composed of uncorrelated regions 
of a size independent of the system size. This can be seen from the following argument: If 
the system is made of many uncorrelated regions, the law of large numbers implies that the 
sum in (|l^) approaches a Gaussian variable of standard deviation proportional to the square 
root of the number of independent zones. Since the size of the zones is independent of L, 
this standard deviation is proportional to L^/^. The average of the modulus of a complex 
Gaussian variable is proportional to its standard deviation, and the normalization factor 
1/A^ = dx/L present in definition (|I8[) completes the factor L'^/"^ in (|19|). This argument 
confirms again the important dynamical role of the 'domains of different wavenumber' iden- 
tified before. To further stablish this picture, we note that the law A ~ L~^/^ should fail 
when there is only one or less than one domain in our system. According to the results of 
Fig. 3c, this will happen for L < 50 (A^ < 200). This has been checked with a system of 
size L = 32 (A^ = 128) so that it is well described by a single local wavenumber during the 
time interval included in Fig. 5b. As seen in that figure, this system does not fit into the 
finite-size-scaling description. A discussion of the dynamics of these small-size systems is 
given in the next section. 

In the dynamics of conventional phase transitions such as order-disorder transitions or 
spinodal decomposition in d > 2 one can also identify independent domains containing some 
'ordered phase'. The difference with our model is that in the usual case such 'ordered phase' 
is close to an equilibrium stable phase. Then domains growth and coalesce until one (or two 
in the case of spinodal decomposition) of them reaches the size of the whole system. This 
domain growth and saturation produces the t/L^ dependence of the scaling function in time. 
In our case, the domains contain an 'unstable periodic phase', so that there is no driving force 
for growth, and the domains keep its size roughly constant (see Fig. 3c), but continuously 
reduce their wavenumber by local roll destruction. This process will presumably continue 
until each domain contains only one of the uniform stable phases ±fe~^/^ and then a growth 



of the domains similar to that in the model (16), as described in |11], is expected 



IV. SMALL SYSTEMS 

Following our discussion above, we will call a system 'small' when its size is smaller than 
the correlation length discussed in the previous section [l^ ~ 50), so that it contains effec- 
tively only one domain of quite uniform wavenumber. In this context it is worth reminding 
here how different characteristics of a system reflect in their Fourier description. First, for a 
small system, modes are more separated than in a larger system {dq = 2n/L). Second, the 
difference between a discrete or a continuous (in x) system is that in the flrst case there is a 
minimum length, dx, so that there appears a maximum wavenumber 2n/dx and only a flnite 
number of modes is at play if L is flnite. In the continuum case {dx — >■ 0), the wavenumber 
cut-off goes to inflnity and we have an inflnity of modes, but their separation continues to be 
determined by L. As usually recognized, a physically continuous system admits a discrete 
description of minimum length dx when 2n/dx is larger than any wavenumber relevant to 
the evolution of the system under study. 

To study the behavior of 'small systems', we have performed numerical integrations of 
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equation (P for a range of values of L for which only three modes are linearly unstable at 
t=0: qn = ndq, n = 0, 1,2. The mode g2 is the most unstable one. The explicit results for 
the amplitude of the three unstable modes and the first linearly stable mode for the case 
L = 18 (A^ = 72, dq ^ 0.349) are shown in Fig. 6a as obtained from the direct numerical 
integration of (|1]) with the initial condition ([T7|). The results indicate that the mode gi 
will dominate the pattern for a long time after an initial short regime in which the fastest 
growing mode g2 dominates. In such small systems, a well defined periodicity is observed 
during a long time interval, so that an apparent nonlinear selection of the mode gi occurs. 
This happens because the stationary unstable solution ipq=q^, dominated by a mode and its 
harmonics, is closely approached during the evolution, in contrast with what we obtain for 
larger systems. Such configuration is, nevertheless, unstable. The pattern finally decays 
and the configuration becomes space homogeneous, as expected. This happens at a time 
t ~ 1.6 • 10^ . This time is not shown in figure 6 and should be compared with the time scale 
in that figure. For such small systems with only three linearly unstable modes, and given 
that the linearly stable modes have a very small amplitude during the whole time evolution, 
it is natural trying to describe the dynamics in terms of three coupled ordinary equations 
for the complex amplitudes ^o, Gqi and 6q^, all the other amplitudes assumed to be zero. 
Truncation to a small number of modes is a largely used technique in the literature OS 



Equation ([l|) is written in Fourier space including only the modes g„ = ndq, with n = 0, 1 
and 2 as 

^,„=o;(g„)e,„-6(a + g„2) ^ O^fi^dq^^q^.q^ (20) 

jj=0,l,2 

The dependence on system size appears through the values of g„, which depend on dq. 

Equation (^) presents stationary solutions whose stability properties are different from 
the corresponding solutions of @). For example, 9q^ = ((1 — gi)/!^^))^''^ and all the other 
amplitudes equal to zero is a stationary solution of (|20D. When the linear stability analysis 
around this solution is carried out, the result depends on the value of dq: The solution is sta- 
ble against homogeneous perturbations if dq < 8~^/^ = 0.354, for any value of a. It is stable 
against perturbations of wavenumber gi only if dq < 7^^/^ = 0.378 (value again independent 
of a). This behavior contrasts with the exact properties of Eq.(p, well reproduced by its 
numerical integration discussed above, for which there are no other stable solutions than the 
one associated to go, independently of the value of dq = 27t/L. The differences in stability 
properties come from the truncation to a small number of modes, equivalent to replacing 
the spatially continuous system by a discretized version. Such differences anticipate that a 
description in terms of a few modes might be qualitatively incorrect. To substantiate this 
point, we have numerically solved Eq. (|20|) with dq = 0.349, the value associated to L = 18, 
and for which the stability properties of the solutions of (|lj) and those of the truncated 



model (^OD are different. The same initial value as in the numerical integration shown in 
Fig. 6a was given to the modes included in (^^). Figure 6b (dotted line) shows that the 
fastest growing mode g2 dominates the final state, in contrast with the numerical integration 
ofEq.(g). 



When Eq.(pOD is enlarged to include the mode gs = 3dq, which is linearly stable at t = 0, 
the stationary solutions in which only one modal amplitude is different from zero are the 
same as above, but its stability properties change. For example, the stationary solution 
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dominated by g2 becomes now unstable for dq > (10/115)2 ^ 0.295 |^. Numerical solution 
of (|20| ) including this fourth mode and keeping dq = 0.349 shows (dashed line in Fig. 6b) 
that 6q2 decays at t ~ 300. Then, a state dominated by gi is approached, as in the integration 
of Eq.(|l|) but, in contrast with that continuous model, this is here the final asymptotic state. 
In this state gs is slightly developed and the other modes have a small amplitude. Again, 
this is a stable solution of the set of 4 equations and the mode g = is not yet reached. 
Figure 6 (solid line) shows that inclusion of up to 7 modes alters the time scales but not the 
quahtative picture. 

As a conclusion, using only the set of linearly unstable modes, or only some additional 
ones, is not enough to describe the time evolution of the continuous system (which we expect 
to be well described by the simulations with 72 grid points, equivalent to 72 modes). A large 
number of linearly stable modes are relevant, although its amplitude remains always very 
small. It is worth noting that sets of equations containing the linearly unstable modes or 
only a few more are often used with success in the literature [|Ii^. The source of success 
in these cases is not usually explicitly stated. It becomes clear from the discussion of the 
example above that the truncation to the set of linearly unstable modes can be useful, at 
most, when these modes are associated with stable stationary solutions of the problem. 
Otherwise, such truncations can produce incorrect results. 
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FIGURES 
FIG. 1. The amplifying factor u>{q) of equation (y) in four different cases. See text for details. 

FIG. 2. Time evolution of 0{x) from a particular initial condition and L = 256. The vertical 
scale is the same for all times, b) Local wavenumber q{x) for the configuration at t = 40 of Fig. 2a. 

FIG. 3. a) Time evolution of the structure factor averaged over 50 runs for L = 128. Times 
showed are, from bottom to top: t = 20, 25, 30, 50, 100, 200, 500 and 750. b) Area A(t) of S{q, t) 
for 4 independent runs (L = 128). c) Correlation length Ic as a function of time (see text for 
details). 

FIG. 4. Number of zeros per unit length in 9{x) versus time for L = 128 (solid line) and 
L = 256 (dashed line). In the inset the same quantity is plotted versus (q) (for L = 128). For 
L = 128, 256 an average over 50, 20 initial conditions respectively was taken. 

FIG. 5. a) Area of the averaged S{q) as a function of time for systems of different sizes. For 
size L = 128 the structure factor has been averaged over 50 different initial conditions. For the 
other system sizes the average is over 20 initial conditions, b) The same area scaled by L^'^. 

FIG. 6. a) Amplitudes 9q^ for n = 0,1,2 and 3 {9q^ being the first linearly stable mode) 
obtained by integrating equation (jl]) for a particular initial condition of the form (|l^ with L = 18. 
b) The same amplitudes obtained from equation (|20| ) (dotted line), from equation ( pO|) enlarged to 
include mode q^ (dashed line), and including up to mode qe (solid line). 
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